Correlation functions of Ising spins on thin graphs 
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We investigate analytically and numerically an Ising spin model with ferromagnetic coupling 
defined on random graphs corresponding to Feynman diagrams of a (j)'^ field theory, which exhibits a 
mean field phase transition. We explicitly calculate the correlation functions both in the symmetric 
and in the broken symmetry phase in the large volume limit. They agree with the results for finite 
size systems obtained from Monte Carlo simulations. 
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I. INTRODUCTION 

The definition of a correlation function on random ge- 
ometries is not so obvious as in the usual fixed geometry 
case. The problem is that we cannot just take two fixed 
points at a given distance because the distance between 
them is changing. One possible solution is to sum over 
all pairs of points at a given distance 



G 



(1) 



Here d{i,j) is the geodesic distance, i.e., the shortest 
path between points i and j. A and B are some quanti- 
ties defined in each vertex whose correlation we want to 
measure. The average is taken over all instances of the 
geometry (configurations) and n is the size of the system. 
Please note that the distance d{i,j) is a very nonlocal 
quantity: in principle it does depend on the whole con- 
figuration not only on the endpoints i and j. Therefore 
the correlation (1) is not a two-point function, which can 
lead to some interesting and non-intuitive behavior. 

Due to this non- locality the calculation of function (1) 
is in general very difficult. As in other fields of statisti- 
cal physics an insight can be gained by examining simple 
solvable models, the best known being probably the Ising 
model. In this paper we study the Ising model on a ran- 
dom geometry ensemble made of all Feynman diagrams 
of a (j)'^ theory. 

The paper is organized as follows: in Sec. II we present 
and analytically solve the model. We proceed on in 
Sec. Ill calculating the spin-spin correlation functions 
using the similarity to Cay ley trees in the infinite vol- 
ume limit. Section IV addresses the influence of spins 
on geometry. Final discussion and summary are given in 
Sec. V. 
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II. THE MODEL 

We consider an Ising ferromagnet spin model on q- 
regular random graphs corresponding to 0* Feynman di- 
agrams of a zero-dimensional field theory [1, 2]. The 
partition function of this model is defined as a sum over 
all 0"^ Feynman diagrams Q and all the values that the 
spins si, . . . , s„ on a graph G can take 







(2) 



Gee si,...,s„ 



If no external magnetic field is present the energy of the 
system on a single diagram G reads 



H{G;si,...,Sn) 



-E ^ 



(3) 



where the sum is over all nearest neighbor pairs. Two 
vertices are considered as nearest neighbors if they are 
connected by a link. This includes loops, in which case a 
vertex is its own neighbor, and multiple links when two 
vertices are counted multiple times in the sum (3). 

Associating with the "up" and "down" spins the 
and 0_ fields respectively, we can generate the requisite 
ensemble from the Feynman diagram expansion of the 
partition function 



Z = d(/)+d0_exp 



1 



where 6^ = 



.) and A is the transfer matrix 



(4) 



(5) 



Following [1] we define the coupling constant g = e^^. 

We will use the saddle-point approximation method to 
calculate the partition function Z in the large n limit. We 
start by performing binomial expansion of the (f)'^ terms 
in Eq. (4) 
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n ,ofc ,q(n— fe) 



n=0 fe=0 ^ ' 



(6) 



Substituting this back into Eq. (4) and taking the term 
for a specific n only we get the partition function of an 
ensemble of graphs with exactly n vertices 

(7) 



where 



C{k,n) 



k\{n-k)V 



(8) 



At this point it is convenient to introduce the rescaled 
fields and the variable z 



k 

z= -, 
n 



and rewrite Eq. (7) as 

Zn=C{n) j^dzC{z,n) j di^+dtj^- e''^'^'l'+''^-\ (9) 

where C(n) is a constant depending on n only and 

, V- ) = - ^ A- V + (? [2 In V+ + (1 - ^) In V-] . 

(10) 

The field integral in Eq. (9) can be asymptotically ap- 
proximated using the saddle-point approximation 



2-K 



n 



(11) 



where J is the Jacobian matrix J^y = /- 



the shorthand notation 



dipxdtpy 

V'± are given by the saddle-point equations 



We use 



(12) 



5Vd 



= 0, 



which can be written in matrix form 



q{l-z)/^_ 



(13) 



(14) 



The above system of two quadratic equations has in gen- 
eral four solutions, from which only the one positive in 
both tp^ and ip- has physical meaning 



ip+iz) = g 



■ 1 + 2^ (5^ - 1) + Vl + 4 {g^ - 1) (1 - z)z 
V^_(^) = V^+(1-^). (15) 
Because of Eq. (14) f{'tl;+,'tl>-) simplifies to 



zlntp+ + {1 — z)il)- 



(16) 



The logarithm of the C{z,n) term can be approximated 
using the Stirling formula 



In (7(2, n) 



So finally 



— n [z In z + (1 — z) In (1 — z)] 

— i In 2(1 — z) + In (27rn) + nlnn — n. 



C{n) / dze"^(^) 



(17) 



(18) 



and combining Eqs. (9), (16), and (17) we obtain F{z) 
to the leading order in n 

F{z) w -^In^- (1 - z)ln(l - z) 

+ q[z\ai)+ + {l- z)\n{l^J\. (19) 

The next non-leading term is given in Appendix A. 

Function F{z) is symmetric around z = 1/2 and in 
the disordered phase it has a maximum at this point. 
Exactly at the transition F"(l/2) = 0, from which we 
can calculate the critical coupling 



9c = 



q-2 



(20) 



For g > gc the F{z) function has a minimum at z = 1/2 
and two maxima at z = 1/2(1 ± m). In saddle-point 
approximation m is the average magnetization 



(M) = 2k-n = n{2z - 1) = nm. 



(21) 



Finding the zeros of the first derivative of F{z) we get 
for g = 3 



5-3 



and 



m 



5-2 V 5 + 1 



(22) 



(23) 



3 



in the g = 4 case. The susceptibihty x = 

can be calculated in saddle-point approximation from the 

integral 



1 



dz i 



leading to 



X 



F"(i(l + m))- 
In the symmetric phase where m — Q wc obtain 

2 9 



X= 1- 



In the broken symmetry phase we get for q — ?i 

4g 



X 



(g -3)(g- 2)2(3 + !)• 



(24) 



(25) 



(26) 



(27) 



Actually, for comparison with Monte Carlo (MC) simu- 
lations we will use 



1 



x^-{{M^)-mY) 



(28) 



instead of x because average magnetization is not 
well defined in numerical simulations. On a finite lattice 
it is in principle zero for all values of /?. However, in the 
broken symmetry phase its measured value will in general 
depend on the algorithm used and the duration of the 
simulation. The average absolute value of magnetization 
is given by the integral 



|A/|) = n 



/p^dz|(2z-l)|e"^(-) 



(29) 



and in the saddle-point approximation in the symmetric 
phase it is equal to 

-{\M\f = -X. (30) 
n TT 

In the broken symmetry phase (|Af|) = nm. We plot 
the resulting expression for x (dashed line) together with 
data obtained from MC simulations in Fig. 1. 

Instead of doing saddle-point approximation of the in- 
tegrals (24) and (29) we can integrate them numerically. 
The results are plotted in Fig. 1 with solid lines. As one 
can see the agreement is very good already for small sys- 
tems even if we keep the leading order terms only given 
by Eq. (19). Including higher order terms does not im- 
prove the result significantly. 



III. CORRELATIONS 

In this section we will calculate the connected spin-spin 
correlation function 




(31) 




FIG. 1. X a function of the coupling g for q — 3. Cir- 
cles (blue) render MC data for graphs of size n = 1000 and 
squares (red) for those ones of size n = 8000 vertices. The 
error bars show statistical uncertainty which increases near 
the transition. The solid lines plot the corresponding results 
of numerical integration whereas the saddle-point approxima- 
tion is given by the dashed line. 



This is not the only way to define the connected correla- 
tion function on random geometry [3]. We will use this 
one because it has the usual property 



X- 



(32) 



Expanding the expression under the average we can ex- 
press (31) as 



= G'''ir)-2mG'^ir)+m^G^\r). (33) 

To calculate the correlation functions we will use the 
fact that in the leading order of n the solution of the 
model coincides with the Bethe solution of the Ising 
model on Cayley trees with coordination number equal 
to q [1]. Contrary to Cayley trees the Ising model on 
Feynman diagrams studied here exhibits a genuine phase 
transition. 

The G^^(r) correlation function is a volume-factor: it 
is the average number of vertices at the distance r from 
a given vertex. On a Cayley tree with fixed geometry it 
is easy to calculate 



G''{r)^q{q-iy 



(34) 
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FIG. 2. Ratio of the correlation function G^^(r) obtained 
from MC simulations of ensembles of various sizes with g = 3 
at p — I3c = iln3 to its infinite volume limit q{q — 1)' ^ . 



In Fig. 2 we plot the ratio of the measured correlation 
function to the infinite volume limit (34). Please note 
the scaling relation evident from the plot 



2G"(r;n) =G"(r + l;2n). 
This can be also written as 



C(r) = nJ" ( — ) = nJ"(rln2 - Inn). 



(35) 



(36) 



We have checked this relation and found that it is very 
well fulfilled already for graphs as small as 1000 vertices. 

The Ising model on Cayley trees is generated by the 
equations 



1 



9-1 



(9-1)! 
1 

(9-1)! VV5 



V9 " 



9-1 



9-1 



(37) 



whose graphical interpretation is shown in Fig. 3. These 
equations are analogous to the ones obtained for Ising 
model on branched polymers [4]. The difference is that 
here the probability of ending the branch is zero. This 
formally means that we have only infinite trees. Hence 
there is no chemical potential associated with every ver- 
tex. 

The partition function can be then obtained as 



1 



Z = - ($^ + $1) 



(38) 



The system (37) will have in general many pairs of so- 
lutions. However, for g < Qc the system is in symmetric 








FIG. 3. Graphical representation of the saddle-point equa- 
tions given by Eqs. (37) for g = 3. The bright bubbles corre- 
spond to "!>+ and the dark ones to Smaller bright and 
dark points stand for vertices carrying the "up" and "down" 
spins respectively. 




FIG. 4. Graphical representation of the correlation function 
Z{r)xy given by Eq. (42) for q — 3. Bubbles depict "I>± and 
points stand for individual vertices connected by links associ- 
ated with A. Only the first and the last vertex have specific 
spins and each of them contributes by a factor. The spin 
of the r—1 vertices along the path connecting them can be 
either "up" or "down" which is covered by the matrix M. 



phase and = The one resulting equation can be 
easily solved yielding 



V9 



1 



(9-1)! 



1 



2 cosh ^ 



(9-1)! 



(39) 



When g > gc the dominant solution of Eqs. (37) will have 
=/= giving a non-zero magnetization [2] 



$9 



(40) 



Finding this solution requires solving the system (37). 
Although this cannot be done analytically for general q 
for 9 = 3 and g = 4 we get in the broken phase 



r V9 



1± 



5-3 



.9 - 1 V V 5 + 1 



for q ~ 3, 



for q = i. 



(41) 



5 



Substituting this into Eq. (40) wc obtain the magnetiza- 
tion which is identical to the expressions (22) and (23) 
calculated using the previous method. 

We will derive the correlation functions by the method 
described in [5]. The correlation function can be repre- 
sented graphically as in Fig. 4. To this picture corre- 
sponds the expression 



Z{r) 



1 



xy 



[(9-1) 



l2 X 



where 



, , _ 1 



(42) 
(43) 



If we rewrite matrix M as 




(44) 



we can see that the correlation functions are related to 
the correlation functions of the Ising spin chain in the 
effective magnetic field {q — 2)/ie// with [6] 



h - 1 , *+ 



(45) 



In the symmetric case in particular heff — 0. 

The required correlation functions can be expressed 
using Z{r)xy as 



G"(r) = ^^Z(r), 



(46) 

G^^(r) = |E2^^W-^' (47) 

xy 

G^'{r) = ^Y.^vZ{rU. (48) 

xy 

After introducing matrices 



( 9-2 

$5^ 



and 



. 1 .... , 1 



(49) 



(50) 



Q = M^AM 

Z{r)xy given by Eq. (42) can be rewritten in the form 



= (9-2)! 



,r y 



{q-l)\ {q-l)V 



(51) 



Assuming that the eigenvalues of Q are Ai and A2, and 
that the corresponding normalized eigenvectors are (a, h) 
and (—6, a) respectively, we can calculate the matrix 
power by diagonalizing Q. 
After some algebra wc obtain: 



1 


AJ(a$| + 


fo*!)^ + A^(6$| -a$i)2 


z 




{q-iy.iq-1) 


1 




-62$l) + A^(62$^-a2$l) 


z 




(9-l)!(9-l) 


1 




6$!)2+A^(6$| + a$i)2 



Z (g_i)!(5_i) 
In Appendix B we show that 
Ai = 9 - 1 , 

(9-1)!^ + " 



(52) 
,(53) 

(54) 
(55) 



A2 = (g - 1) 
= (9-1) 



(9-1)! 
(9 - 1) A2 , 



)-l 

cosh((g - 2)heff) - 1 



where 



The corresponding eigenvector is equal to 
(a, 6) = , 

Substituting this into Eq. (52) we get 

1 {q-iy-^ 



G-(r) 



Z (g-1)! 



(56) 
(57) 

(58) 

(59) 



Using the definition (38) of Z we finally obtain the result 
(34) which is a check of the consistency of the method 
used. 

For spin-spin correlation functions we obtain 



G^^r) = 



G''{r) 



^, {q - 1)' ($; - $1)2 + 4$«_$1A^ 

T\2 



(g-l)!(g-l)($^ + $«) 



G"(r) 
G"(r) 



4$^$1 

($9^ + $9)2 

cosh {q heff) 



A5 -hm^ 



(61) 



In the symmetric case when m = and = are 
given by Eq. (39) we readily get 



5-(r) = 



G-'^r) 
gTi(0 



.9-1 
5+1 



tanh''/?. (62) 
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FIG. 5. Correlation function g'"'{r) for q = 3. Symbols render 
MC data for graphs of size n = 256000 in the symmetric phase 
for P — 0.45 (circles) , at the transition /3 = /3c = | In 3 (squares) , 
and in the broken phase for /3 = 0.60 (diamonds). The solid 
lines plot the analytical predictions given by Eqs. (62) and 
(63). 



This is as predicted the result obtained for correlation of 
Ising spins on the chain [6]. It is easy to check that the 
relation (32) is satisfied by the above function. 

In the broken phase inserting Eqs. (41) into Eq. (60) 
we obtain for 5 = 3 



1 



and for q = 4 



(g-2)2(g + l) (g-iy 



(52 - 2)2 (52 _ 1). 



(63) 



(64) 



Although similar formulas can be derived for higher val- 
ues of q they are much more complicated. 

Figure 5 plots the correlation function (r) for g = 3 
and different values of /3. As one can see the agreement 
with MC results is very good. In Fig. 6 correlation func- 
tions for different g's at the respective transition points 
I3c{q) = ^ In compared. Again, MC results match 

the asymptotic predictions. The small discrepancies di- 
minish with the increase of the graphs' size. 



IV. INFLUENCE OF SPINS ON GEOMETRY 

On a finite tree spins can be integrated out exactly and 
the resulting factor does not depend on the shape of the 
tree. This means that the spins do not have any influence 
on the geometry. The situation changes when cycles are 
allowed, which is especially obvious for loops, i.e., links 



FIG. 6. Correlation function g'^ir) at the transition I3 = l3c{q) 
for different values of q. Points represent MC results for 
graphs of size n — 64000 while the solid lines plot the ana- 
lytical predictions. Please note that for q — 6 the maximum 
encountered distance (the diameter of the graph) is r = 9. 
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FIG. 7. Average number of loops n; in graphs of size n = 32000 
vertices and g = 3. The black points represent MC data and 
the dashed red line marks the transition point /? = /Re- 



attached at both ends to the same vertex. Without any 
spins or /3 = one can show combinatorially that the ex- 
pected number of such loops on a graph equals {q—l)/2. 
For positive /3, however, we should observe an enhance- 
ment as each loop contributes the factor to the par- 
tition function contrary to links joining different vertices 
which can have different spins. For /? — >■ 00 all the spins 
have the same sign and again the geometry decouples. 
In Fig. 7 we plot the average number of loops ni as a 



7 



function of /3. The results agree qualitatively with the 
above scenario. Nevertheless, one should note that while 
looking pronounced this is still only a 1/n effect. The 
number of loops is independent of n and negligible in the 
large volume limit. 



V. SUMMARY AND DISCUSSION 

We have analyzed in detail a simple model of spins on 
a random geometry. In the large n limit it is formally 
equivalent to the Ising model on an infinite Cayley tree. 
However, it is well defined, has a genuine phase tran- 
sition, and can be easily simulated using Monte Carlo 
methods. 

We have derived expressions for correlation functions 
in both the symmetric and the broken phase with meth- 
ods less formal then in [6]. From these calculations we 
obtain a picture of the transition. The correlation func- 
tion g^^{r) does not exhibit any critical behavior: the 
correlation length is finite for any finite /3. Neverthe- 
less, the volume factor G^^(r) grows exponentially and 
this growth offsets the decay of the correlation function 
g^^{r). As we increase r by one, the influence of the spins 
at this distance drops by {g — l)/{g + 1), however, the 
number of these spins increases by g' — 1. When 



5+1 



(65) 



all shells of spins contribute equally and we observe the 
phase transition. It is easy to check that this gives the 
expression (20) for the critical value of g. 
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Appendix A: Next to leading corrections to F{z) 
The Jacobian matrix J^y = f-:;p :^ is equal to 

'z/^+ 



J = A-i 



(1 - z)/i;_ 



(Al) 



from which 



det J 



\V'+ V'- y 



o z(l — z) 

+1 J - '2 - 



(^-V'+; 
(A2) 

and the approximation of the function F{Z) given by 
Eq. (19) takes the form 

F{z) w - [z\nz + {l - z) ln(l - z)] 
+ q[z ln'0+ + (1 — 2;) In-i/J^] 
- 5^ {In (det J) + In [^(1-^)]}. (A3) 

Appendix B: Eigenvalues 



Equations (37) can be rewritten as 
1 /^/5^r' 



V9- 



$4 



(g-1)! [j^^l-^ V5*rV V^- 



Q 



(Bl) 



The above equation means that matrix Q has an eigen- 
value equal to one and that the corresponding eigenvector 
is proportional to ($+,$_). 

Matrix Q defined in Eq. (50) can be expressed as 



Qxy — (<? 1) Qxy ( ^ 



(B2) 



One can check that this implies that Q has an eigenvalue 
Al = g' — 1 and the corresponding normalized eigenvector 



(a, 6) = 



(B3) 



The second eigenvector is perpendicular to this one. Mul- 
tiplying it by Q we obtain the second eigenvalue. 
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